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Abstract 

We report branching fractions and CP asymmetries for B — > nK* and B — ► np decays. These 
results are obtained from a 414 fb -1 data sample collected at the T(4S) resonance, with the Belle 
detector at the KEKB asymmetric energy e + e~ collider. The branching fractions, in parts per 
million, of nK*°, nK* + , np°, and np + are 15.9 ± 1.2 ± 0.9, 19.7±f;g ± 1.4, 0.84±g;ff ± 0.18, and 
4.ll 1 ; 3 ± 0.34, respectively. We find no evidence for CP asymmetries in these modes. 

PACS numbers: 13.25.HW,14.40.ND 
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INTRODUCTION 



Charmless hadronic B decays play an important role in understanding CP violation in 
the B meson system. Studies of B — > rjK* and B — > r/p lj are important examples of such 
decays. In the standard model (SM), Penguin (tree) diagrams are expected to dominate in 
B — > rjK* (B — > rjp) decays. The large branching fraction for B — > rjK* compared to that 
for for B -> rjK decay 0, S, 3 can be explained qualitatively in terms of the interference 
between non-strange and strange components of the r\ meson, but are higher than recent 
theoretical predictions [H|, E| 0, El • In a similar vein, the larger measured branching fraction 
for charged (B + — > 7]K* + ) vs. neutral (B° — > f]K*°) decays may suggest an additional 
SU(3)-singlet contribution 0, 0, E| or constructive interference between SM penguin and 
tree amplitudes or between SM- and New Physics(NP)-penguin amplitudes. 

In the standard model, direct CP violation (DCPV) occurs in decays that involve two (or 
more) amplitudes that have different CP conserving and CP violating phases. The partial 
rate asymmetry can be written as 
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where / denotes a self-tagging final state, B is either B + or B° meson, B and / are the 
conjugate states, and AS (A0) is the difference of the CP-conserving (CP- violating) phases 
between amplitudes A\ and A 2 . The charge asymmetry will be sizeable when the two ampli- 
tudes are of comparable strength with significant phase differences. SM Penguin diagrams 
are expected to dominate in the B — > rjK* decays, while tree diagrams are expected to 
dominate in B + — > i]p + decays. Therefore, DCPV in both cases should be small. However, 
SM- and NP-penguin diagrams may interfere to generate sizable direct CP violation. 



DATA SET AND APPARATUS 

This analysis is based on a data sample collected at the T(4S) resonance with the Belle 
detector [13] at the KEKB accelerator. The data sample corresponds to an integrated 
luminosity of 414fb _1 and contains 449 x 10 6 BB pairs. 

The Belle detector is designed to measure charged particles and photons with high ef- 
ficiency and precision. Charged particle tracking is provided by a silicon vertex detector 
(SVD) and a central drift chamber (CDC) that surround the interaction region. The charged 
particle acceptance covers the laboratory polar angle between 9 = 17° and 150°, measured 
from the z axis that is aligned anti-parallel to the positron beam. Charged hadrons are 
distinguished by combining the responses from an array of silica aerogel Cherenkov counters 
(ACC), a barrel-like array of 128 time-of-fright scintillation counters (TOF), and dE/dx 
measurements in the CDC. The combined response provides K/tt separation of at least 
2.5cr for laboratory momentum up to 3.5 GeV/c. Electromagnetic showers are detected in 
an array of 8736 CsI(Tl) crystals (ECL) located inside the magnetic volume, which cover 
the same solid angle as the charged particle tracking system. The 1.5-T magnetic field is 
is contained via a flux return return that consists of 4.7 cm thick steel plates, interleaved 
with resistive plate counters used for tracking muons. Two inner detector configurations are 
used. A 2.0 cm beampipe and a 3-layer silicon vertex detector are used for the first sample 
of 152 x 10 6 BB pairs, while a 1.5 cm beampipe, a 4-layer silicon detector and a small-cell 
inner drift chamber are used to record the remaining 297 x 10 6 BB pairs [l2T ]. 
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For Monte Carlo (MC) simulation study, the signal events, generic b — > c decays and 
charmless rare B decays are generated with the EVTGEN k| event generator. The contin- 
uum MC events are generated with the e + e~ — > 7* — > qq process in the JETSET generator. 
The GEANT3 [13] package is used for detector simulation. 



EVENT SELECTION AND RECONSTRUCTION 

Hadronic events are selected based on the charged track multiplicity and total visible 
energy sum, which give an efficiency greater than 99% for BB events. All primary charged 
tracks are required to satisfy track quality cuts based on their impact parameters relative 
to the run-dependent interaction point (IP): within ±2 cm along the z axis and within 
±1.5 cm in the transverse plane. Particle identification (PID) is based on the likelihoods 
£k/(At + Lk) for charged kaons and pions, respectively. A higher value of £k/(At + Lk) 
indicates a more kaon-like particle. PID cuts are applied to all charged particles except 
pions from K s in this analysis. Unless explicitly specified, the PID cuts are £k/(At + Lk) 
> 0.6 for kaons and < 0.4 for pions. The PID efficiencies are 85% for kaons and 89% 
for pions, while the fake rates are 8% for pions faking kaons and 11% for kaons faking 
pions. In forming tt° candidates from photon pairs, the photon energies must exceed 50 
MeV and the 7r° momentum in the center of mass (CM) frame must exceed 0.35 GeV/c. 
Kg candidates are reconstructed from pairs of oppositely charged tracks whose invariant 
mass lies within ±10MeV/c 2 of the K s meson mass. We also require the vertex of the K s 
to be well reconstructed and displaced from the interaction point, and the K s momentum 
direction be consistent with the K s flight direction. 



i] Meson Reconstruction 

Candidate 77 mesons are reconstructed through 77 — > 77 and 77 — > 7r + 7r~7r°. If one of the 
photons from the former 77 decay mode could be paired with another photon with recon- 
structed 77 mass within 3a of the 7r° meson mass, then the 77 candidate is discarded. We 
relax the PID requirement for charged pions from the latter 77 decay mode to Ck/ (At + -^k) 
< 0.9. The momenta of the 77 candidate is recalculated by applying the 77 mass constraint for 
B reconstruction. The 77 — > 77 candidates must satisfy | COS6 1 *! < 0.90, where 9* is the angle 
between the photon direction in the 77 rest frame and the 77 momentum in the CM frame, 
so as to suppress the soft photon combinatorial background and B — > K*^ feed-across. The 
reconstructed mass resolutions are 12 MeV/c 2 for 77 — ► 77 and 3.5 MeV/c 2 for 77 — > 7r + 7r~7r°. 



K* , p Meson Reconstruction 

K*° candidates are reconstructed from K~ir + and A^7r° pairs, while K* + mesons are 
reconstructed from K + tt° and K s tt + pairs. Candidate K* mesons are required to have re- 
constructed masses within ±75 MeV/c 2 of the nominal value. Candidate p° (p + ) mesons are 
reconstructed from tt^tt + (tt ti + ) pairs. Each combination is required to have a reconstructed 
mass within ±150 MeV/c 2 of the nominal value. 
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B Meson Reconstruction 

B meson candidates are reconstructed from t]K*°, t]K* + , rjp , and rjp + combinations. 
They are characterized by the beam- constrained mass M bc = \J E^ cam — \Pb\ 2 and the energy 
difference AE = E B — E hcam , where E heam = 5.29 GeV, and Pb and E B are the momentum 
and energy, respectively, of B candidate in the T(4S) rest frame. We define the fit region in 
the M hc -AE plane as M bc > 5.2 GeV/c 2 and \AE\ < 0.25 GeV. We define the signal region 
as the overlap of the bands M bc > 5.27GeV/c 2 and \AE\ < 0.05 GeV. 

One candidate per event is required for all modes. The best candidate is chosen based on 
the sum of r\ vertex mass constraint fit x 2 an d K*°(p°) vertex x 2 for B — > 7]K*°, K*° — > K + tt~ 
and B — > r)p°. For the B + — > r]K* + , K* + — > K°7T + mode, the reconstructed K* + mass x 2 
and the r\ vertex mass constraint fit x 2 are used. For the remaining modes, the x 2 values of 
the 7] and n° mass constraint fits are used. 

BACKGROUND SUPPRESSION 

The dominant background for exclusive two-body B decays comes from the e + e _ — > 
7* — > qq continuum (q = u, d, s, c), which has a jet-like event topology compared to the 
spherical BB events. Other major backgrounds involve feed-across from other charmless 
B decays. For decays involving the p meson, an additional background arises from b — >• c 
decays. However, the impact of this background is small since the M^ c and AE distributions 
do not peak in the signal region. 

Continuum Background 

Signal and continuum events are distinguished by two means. First, we require | cos^tI < 
0.9, where 9t is defined as the angle between the rj daughter of a B candidate and the thrust 
axis from all the particles in the event not associated with that B candidate. This retains 
90% of signal and removes ~56% of continuum. Second, a probability density function C s 
(£ c ) for signal (continuum) is formed from two independent variables — cos#g, where 8b is 
the polar angle of the B candidate, and a Fisher discriminant T = a ■ R that combines 
seven event shape variables: cos 6t, S± (the scalar sum of the transverse momenta of all 
particles outside a 45° cone around the B candidate direction divided by the scalar sum of the 
momenta), and the five modified Fox- Wolfram moments Rf, R%°, R^°, and R™. The 

Fisher discriminant's weight vector a is determined by optimizing the separation between 
signal events and continuum background using MC data; these Fox- Wolfram moments are 
used since they are not correlated with M bc . The likelihood ratio CTZ = C S /(C S + C c ) is 
used to distinguish signal from continuum, since it peaks near 1 for signal and near for 
continuum. 

The distribution of CTZ is found to depend somewhat on the event's B flavor tagging 
quality parameter r (l7j . which ranges from zero for no flavor identification to unity for 
unambiguous flavor assignment. We partition the data into three r regions, and determine 
the optimal cut on CTZ in each region by maximizing the significance Ng / \J N$ + Nb formed 
from the retained number of signal (Ns) and continuum background (Nb) events in MC 
samples. A typical cut of CR > 0.5 is ~ 75% efficient for signal and removes ~ 81% of the 
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continuum background, while CR > 0.9 is ~ 38% efficient for signal and removes ~ 97% of 
the continuum background. 

Feed-across and b — > c Backgrounds 

A B — > K*(p)"f veto is applied to suppress the K*(p)^ feed-down for B — > r)K*(p), 
7] — > 77 decays. The MC studies show that the background from b — *> c decays and the 
feed-across from other charmless B decays are negligible for B — > 77/^* decays. For B — > r^p 
decays, the background can be as large as 4% of the total yield (Table H}. The contributions 
of these backgrounds are taken into account in the analysis. 

TABLE I: Estimated yields from b — > c (iVb c ), charmless B decay (iV r ) backgrounds, r/X* feed- 
across (A^feed) and actual yields from all sources (N) after the cos 9t and CTZ cuts. 



Mode 


N hc N T 




N 




62 81 


17 


2931 




67 27 


5 


1063 




148 74 


3 


4169 


VnTCTT^ P 


76 22 


1 


1809 



ANALYSIS PROCEDURE 

Signal yields are obtained using an extended unbinned maximum likelihood (2-D ML) fit 
to the Mbc and AE distributions in the M\> C -AE fit region that satisfy the cos 9t and CTZ 
requirements. 

For iV input candidates, the likelihood is defined as 

(Ns+NB+Nbc+Nl+N{ccd) N 

L(N S , N B ) = Y[[N s P Si + N B P Bi + N hc P hCi + N v P ri + N ieed P {eedi ) (2) 

where Ps 4 , P Bi , Pb Ci , and Pf ee di are the probability densities for event i to be the signal, 
continuum, b — *> c, charmless B decay and feed-across backgrounds for variables M^ c and 
AE, respectively. Poisson statistics for Ns and N B , the extracted yields for the signal and 
continuum background from the fit, are considered in this type of likelihood L. The yields 
Abe, A r , and Af ee d are fixed from the MC analysis. 

The continuum, b —* c and charmless B-decay background AE PDF's are modeled by 
higher order polynomial functions. The continuum and b — > c background components in 
Mb c are modeled by a smooth function Due to the peaking behavior of M^ c in the 

signal region from charmless B decay backgrounds, we use a sum of two bifurcated-Gaussian 
functions to model the distributions. The bifurcated Gaussian combines the left half of a 
wide-resolution Gaussian with the right half of a narrow-resolution Gaussian, both having 
a common mean. For B —>■ rjp decays, M^ c and AE distributions from 7]K* feed-across 
will behave like signal with a AE shift of —50 MeV. The PDF shape for each contributions 
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is determined by MC. The slope of the continuum-background AE polynomial and the 
parameters of the Mt, c function are allowed to float in each fit. 

For signal AE, we used two bifurcated-Gaussian functions. The first accounts for 60- 
80% of the total area and the second, which has a larger width, is used to model the 
low-energy tail. M^ c is weakly correlated with AE, so we construct separate bifurcated 
Gaussians for M bc in the three ranges \AE\ < 0.05 GeV, 0.05 GeV < \AE\ < 0.1 GeV, and 
0.1 GeV < \ AE\ < 0.25 GeV. 

For decays with more than one sub-decay process, the final results are obtained by fitting 
the sub-decay modes simultaneously with the expected efficiencies folded in and with the 
branching fraction as the common output. This is equivalent to summing \ 2 — — 21n(L) as 
a function of the branching fraction for the sub-decay channels. The statistical significance 



(£) of the signal is defined as y — 2 \n(L / L maiX ) , where L and L max denote the likelihood 
values at zero signal events and the best fit numbers, respectively. 

The 90% confidence level (C.L.) upper limit x 90 is calculated from the the equation: 

C#T^r = 90% • ( 3 ) 

For this calculation of X90, the likelihood function is modified to incorporate the systematic 
uncertainty. 




FIG. 1: Projections on M^ c (for the signal slice in AE) and AE (for the signal slice in M^, c ) for 
rjK*° (a,b) and rjK*^ (c,d) with the expected signal and background function overlaid. The shaded 
area represents rj — > tt + tt^-k° decays. 



MEASUREMENTS OF BRANCHING FRACTIONS 

The overall reconstruction efficiency e is first obtained using MC samples and then mul- 
tiplied by PID efficiency corrections obtained from data. The PID efficiency correction is 
determined by using D* + — ► D°7T + , D° — > K~n + data samples. 
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TABLE II: Summary of results for each channel listed in the first column. The measured signal 
yield (Ns), statistical significance (£), reconstruction efficiency (e), total efficiency (etot) including 
the secondary branching fraction, and measured branching fractions are shown. Uncertainties 
shown in second and sixth columns are statistical only. The upper limit for r]p° is at 90% C.L. 
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Other MC efficiency corrections are determined by comparing data and MC predictions 
of other well-known processes. The charged tracking efficiency correction is studied using 
a high-momentum 77 sample, comparing the ratio of 77 — > 7r + 7r~7r° to 77 — > 77 between data 
and MC. The same high-momentum 77 sample is also used for rr reconstruction efficiency 
corrections by comparing the ratio of 77 — > 7r 7r 7r° to 77 — > 77 between the data and MC sam- 
ple. The Kg reconstruction efficiency is verified by comparing four K*(892) decay channels 
(K + ir~, K + tt°, K s tt + , K S 7T°) in inclusive K* and exclusive B — ► J/ipK* samples. The £i? 
cut efficiency correction is determined using B — > Dir + decays. For 77 and K* reconstruction 
and mass cuts, we use the high-momentum 77 and K* sample for the efficiency correction 
studies. 

The above studies show good agreement between the efficiencies in the data and MC 
sample at about the 2% level. The PID, rr , 7/ and K* reconstruction efficiency corrections 
are applied and the systematic uncertainties are also obtained from the above studies. 

The fitted signal yields and branching fractions are shown in Table ILT1 Several consistency 
checks are made, including tighter CR cuts, 1-D ML M^ c , and AE fits, and they are all shown 
to be consistent. The total observed yields from the fits, with all sub-decay modes combined, 
are N nK *o = 459.2±|f;l for B° -> r]K*°, N vK ,+ = 253.4l|^ for B + -> r]K* + , N vp o = 20.4tg;jj 
for B° -> rip and N vp + = 53.91 \f f for B + — > 7]p + . Figure [U shows the projections of of the 
data and the fits onto M^ c (for events in AE signal slice) and AE (for events in the M^ c 
signal slice) for the B — > 77^* decays, while Fig. |21 shows the corresponding projections for 
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FIG. 2: Projections on M\, c and AE from 2-D ML fitting results for 7jp° (a,b) and rjp + (c,d) with 
the expected signal and background function overlaid. The shaded area represents rj — > 7r + 7r~7r° 
decays. 

the B rjp decays. For 5 + — > r^p + decays, where a clear excess is seen, we examine the 
properties of the p + candidates. Although statistically limited, clear p + mass peaks and a 
polarized cos #hei distribution are observed (Fig. |3J) and are consistent with the expectation 
from r]p + with no significant non-resonant B + — > r]7r + 7r . 



SYSTEMATIC ERROR 

Systematic errors arise from efficiency corrections and fitting. The main sources of uncer- 
tainties in the efficiency corrections are from the reconstruction of low-momentum charged 
tracks, low-energy photon finding, and CR cut efficiency, each at a level of a few per cent. 
The systematic errors include contributions of 1% for CR cuts, 1% per reconstructed charged 
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FIG. 3: Fitted yields vs (a) p + helicity, (b)7r + 7r° invariance mass from B — > 7]p + decays. The 
overlaid histograms are expected distributions from MC and normalized by the 2-D fit results. 
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FIG. 4: Fitted yields vs the Kir invariant mass for the (a) K*° and (b) K* + modes. Overlaid 
functions are results from fitting with a K* line shape and a higher resonance. Only the parameters 
of the higher resonance are free during the fit. 



particle, 0.5% for each charged particle identification, 4% for 7r° reconstruction, 4.5% for Kg 
reconstruction, and 2% for r\ reconstruction with t] — > 77. To evaluate a possible non- 
resonant or higher resonance contribution in the K* region, we perform a 2D ML fit within 
different Ktt invariant mass region from 0.7 GeV/c 2 to 1.72 GeV/c 2 . Figure 0] shows our 
results with a K* line shape and a higher resonance overlaid with the K* area fixed. Based 
on this study, the non-resonant Kn contributions are 1.7 ± 0.5% for t]K*° and 3.7 ± 2.0% 
for i]K* + decays. These corrections will be applied to the final branching fraction measure- 
ments of B — > r/K*. Due to the limited statistics for B + — > i]p + decays, a larger systematic 
error for the non-resonant or higher resonance contributions is assigned with no corrections 
applied. 

We use B° — > D°7r + decays to estimate the uncertainties in the signal function PDF's 
used for fitting Mt, c and AE by comparing the mean and the width of the M^ c and AE 
distributions between the B° — > D°tt + data and the MC sample. For MC estimated b — > c 
and charmless B decay backgrounds, we vary the estimated yields by 50% and refit the data. 
The difference in the fit yields from the nominal values contributes to the systematic error. 
The overall relative systematic errors are 5.8% for t]K*°,7.2% for r)K* + , 22% for r/p and 
8.4% for 7]p + . 



A C p MEASUREMENTS 

We measure Acp in the decays of B — > rjK* and B + — > i]p + . Due to the wrong-tag 
fraction w, the value of Acp is not the same as the measured values A^p. Their correlation 
is Aq P = (1 — 2w)Acp- In the decay modes we study, only those in which Aqp values are 
determined by low momentum charged pions will have a non- negligible w. The wrong-tag 
fractions for K* + — > K°tc + is ~ 1.5% for t]K* + decays and ~ 2.0% for r)p + decays. Other 
decays have w < 0.1%. Since the result for 7]K* + is obtained from a simultaneous fit to all 
four sub-decay modes with roughly equal statistics for K* + —>■ K + tt° and K* + — > K°n + , 
the wrong-tag effect for Acp should be < 0.7%. Thus, there are no wrong-tag fraction 
corrections for 7]K* decays. For i]p + , 2% is assigned for w. 

To incorporate the CP asymmetry in the fit, the coefficients of the signal and continuum 
background PDFs in the likelihood are modified: N s — > |AT 5 (1 — qA^p) and Np — > ^Nb(1 — 
qAcp,qq), where q = +1(— 1) for a B(B) meson tag and A°p S ,A cpm are the A cp outputs for 
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FIG. 5: Projection on Mb c (left) and AE'(right) from 2-D ML fitting results for r]K* with the 
expected signal and background function overlaid. Figures (a,b,e,f) are for b decays and (c,d,g,h) 
are for b decays. 

signal and continuum respectively. The results are A°p S (r]K*°) = 0.17 ±0.08, A°p S (r]K* + ) = 
0.03 ± 0.10 and A°p S (rjp + ) = — O.OCtoJa- Figure El shows the projections of of the data and 
the fits onto M^ c (for events in AE signal slice) and AE (for events in the M^ c signal slice) 
for the B — > r]K* decays. 

Since the systematic errors in r\ reconstruction and the number of BB events cancel in the 
ratio, the systematic uncertainty of Aqv is mainly from the asymmetry for PID of charged 
kaons and the fitting PDF's. The efficiency asymmetry for the PID of charged kaons is 0.01 
in absolute value. 



DISCUSSION AND CONCLUSION 

In summary, we report measurements of the exclusive two-body charmless hadronic B — > 
rjK* and B — > rjp decays with high statistics. Our results confirm that the branching 
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fractions for B° — > r]K*° and B + — > rjK* + are large and are consistent with previous 
measurements @, 0|- The branching fractions obtained are B[B — > r/K* ) = (15.9 ± 1.2 ± 
0.9) x 1(T 6 , and B(B -> r/K* + ) = (19.7i?;g±1.4) x 10~ 6 , where the first error is statistical and 
the second systematic. Our measurements imply that the branching fraction of B + — > i]K* + 
is la higher than B° — > f]K*°, which may constrain the contribution of the flavor singlet 
penguin amplitude as suggested by various theorists [1, H|. Large signals of £> + — > r/p + 
have been observed with the branching fraction B{B — > r]p + )= (4.ll}; 3 ± 0.34) x 10~ 6 . The 
branching fractions and 90% C.L. for B° -> r?p decays are B(B° r/p°)= (0.84±g;|f ± 
0.18) x 10" 6 (< 1.9 x 10~ 6 ). 

We also search for direct CP asymmetry in B — > r^fC and 5 + — > i]p + , and final results are 
consistent with no asymmetry. The results are: A cp {i]K* ) = 0.17±0.08±0.01, A C p{f]K* + ) = 
0.03 ± 0.10 ± 0.01, and A cp {r]p + ) = -0.04±Jj;H ± 0.01. 
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